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, Abstract. We construct different integrable generalizations of the massive Thirring cqua- 

' tions corresponding loop algebras g"^ in different gradings and associated "triangular" R- 

, operators. We consider the most interesting cases connected with the Coxeter automor- 

^ ■ phisms, second order automorphisms and with "Kostant-Adler-Symes" i?-operators. We 

recover a known matrix generalization of the complex Thirring equations as a partial case 
of our construction. 
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1 Introduction 

> \ 

'st" \ A theory of hierarchies of integrable equations in partial derivatives is based on the possibility 

to represent each of the equations of the hierarchy in the so-called zero-curvature form: 



dU(x,t,u) dV(x,t,u) , , , 

(N; ^^ + [C/(x,t,n),F(x,t,n)] =0, (1) 

00 

\ where U (x, t, u), V{x, t, u) are g- valued functions with dynamical variable coefficients, q is simple 

(reductive) Lie algebra and u is an additional complex parameter usually called spectral. In order 
for the equation ([1]) to be consistent it is necessary that U{x, t, u) and V{x, t, u) belong to some 
\ closed infinite-dimensional Lie algebra g of g-valued functions of u. 

d ' There are several approaches to construction of zero-curvature equations ([1]) starting from 

Lie algebras g. All of them are based on the so-called Kostant-Adler-Symes scheme. One of 
the most simple and general approaches is the approach of [3l [10] and [H [7] that interprets 
equation ([T]) as a consistency condition for a two commuting Hamiltonian flows written in the 
Euler-Arnold or Lax form. In the framework of this approach elements U{x, t, u) and V{x, t, u) 
coincide with the algebra- valued gradients of the commuting Hamiltonians constructed with the 
help of the Kostant-Adler-Symes scheme, the cornerstone of which is a decomposition of the 
Lie algebra g (as a vector space) in a direct sum of its two subalgebras: 

5 = 0+ + 0-- 

The algebra-valued gradients of the commuting Hamiltonians coincide with the restrictions of 
the algebra- valued gradients of Casimir functions of g onto the subalgebras q±. Hence, such 
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approach permits [H] to construct using Lie algebra q the three types of integrable equations: 
two types of equations with U-V pair belonging to the same Lie subalgebras q± and the third 
type of equations with [/-operator belonging to 0+ and ^-operator belonging to g_ (or vise 
verse). The latter equations are sometimes called "negative flows" of integrable hierarchies. 

Nevertheless the approach of [3l [U [7] does not cover all integrable equations. In particular, 
it was not possible to produce by means of this approach integrable equations possessing U-V 
pairs with [/-operator belonging to 0+ and F-operator belonging to 0_ in the case 0+ n0_ 7^ 0, 
i.e. in cases dropping out of the scope of Kostant-Adler-Symes method. An example of such 
situation is the Thirring integrable equations with the standard s/(2)-valued U-V pair [9]. 

In the present paper we generalize the method of [3l [10] of producing soliton equations and 
their U-V pairs, filling the gap described above, i.e. making the method to include all known 
soliton equations, and among them Thirring equation. We use the same idea as in [3] and [6], 
i.e. we interpret zero-curvature condition as a compatibility condition for a set of commuting 
Hamiltonian flows on 0* but constructed not with the help of Kostant-Adler-Symes method 
but with the help of its generalization - method of the classical i2-operator, where R : q ^ g, 
satisfies a modified Yang-Baxter equation [12]. In this case one can also define [12] two Lie 
subalgebras qr^ such that qji_^ + 0/j = 0, but in this case this sum is not a direct sum of two 
vector spaces, i.e. qr^ n Qfi_ ^ 0. 

Hence, in order to achieve our goal it is necessary to construct with the help of the R- 
operator an algebra of mutually commuting functions. Contrary to the classical approach of |12j 
they should commute not with respect to the ii-bracket { , }/j but with respect to the initial 
Lie-Poisson bracket { , } on 0*. In our previous paper [13] the corresponding functions were 
constructed using the ring of Casimir functions I'^{q*). In more detail, we proved that the func- 
tions I^^{L) = Ik{{R* ± 1){L)), where Ik{L),Ii{L) e I^(q*), constitute an Abelian subalgebra 
in C°°{q*) with respect to the standard Lie-Poisson brackets { , } on 0*. The algebra- valued 
gradients of functions the I^"^ (L) belong to the subalgebras correspondingly. 

In the case when the i2-operator is of Kostant-Adler-Symes type, i.e. R = — P-, where 
P± are projection operators onto subalgebras qr^ = g± and 0+n0_ = we re-obtain the results 
of [To] (see also [13]) as a partial case of our construction. In the cases of more complicated 
iJ-operators our scheme is new and generalizes the approach of [10]. In particular, the important 
class of U-V pairs satisfying zero-curvature equations that can be obtained by our method are 
connected with the so-called "triangular" i?-operators. In more detail, if possess a "triangular" 
decomposition: = 0+ + 0o + 0-, where the sum is a direct sum of vector spaces, q± and 0o 
are closed Lie subalgebras and g± are 0o-modules, Rq is a solution of the modified Yang-Baxter 
equation on 0o, P± are the projection operators onto the subalgebras q± then R = P+ + Rq — P^ 
is a solution of the modified Yang-Baxter equation on (see [13] for the detailed proof). The 
Lie subalgebras qr^ have in this case the following form: gRj_ = 0± + 0o> i-e. Qr_^ H gR_ ^ 0. 
Such i?-operators are connected with the Thirring-type integrable models. These i?-operators 
were first considered in [1] and in [5], where the usual sZ(2)-Thirring equation was obtained 
using "geometric" technique. 

In order to obtain the Thirring integrable equation and its various generalizations in the 
framework of our algebraic approach we consider the case when coincides with a loop algebra 
= (B) Pol{u~^ , u) or its "twisted" with the help of finite-order automorphism a subalgeb- 
ra 0°^. The algebras 0°^ possess the natural "triangular" decomposition with the algebra 0o being 
a reductive subalgebra of stable under the action of a. For each algebra 0°^ with a natural 
triangular decomposition and for each classical R operator Rq on 0o we define an integrable 
equation of the hyperbolic type which we call the "non- Abelian generalized Thirring equation" . 
We show that in the case = s^(2), = 1 and Qo = i), where f} is a Cartan subalgebra of s/(2) 
our construction yields the usual Thirring equation and its standard sZ(2)-valued U-V pair. 
We consider in some detail the cases of the generalized Thirring equations that correspond to 
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the second order automorphism of Q. For the case of such automorphisms and g = gl{n) and 
special choice of Rq we obtain the so-called matrix generalization of complex Thirring equations 
obtained by other method in [15]. After a reduction to real subalgebra u{n) these equations 
read as follows: 



where ^± G Mat(p, q), n± G M are constants, i.e. are exact matrix analogues of the usual massive 
Thirring equations. 

We also consider in detail the case of the generalized Thirring equations that correspond to 
the Coxeter automorphisms of q. We call the generalized Thirring equations corresponding to 
the the Coxeter automorphisms the "generalized Abelian Thirring equations". We show that 
the number of independent fields in the generalized Abelian Thirring equation corresponding to 
Coxeter automorphism is equal to 2(dimg — rankg) and the order of the equations grows with 
the growth of the rank of Lie algebra g. We consider in detail the generalized Abelian Thirring 
equations corresponding to the case g = sl{3). 

For the sake of completeness we also consider non-linear differential equations of hyperbolic 
type corresponding to and the Kostant-Adler-Symes /^-operator. We show that the obtained 
equations are in a sence intermediate between the generalized Thirring and non-Abelian Toda 
equations. The cases of the second order and Coxeter automorphisms are considered. The sl{2) 
and sl{3) examples are worked out in detail. 

The structure of the present paper is the following: in the second section we describe com- 
mutative subalgebras associated with the classical R operator. In the third section we obtain 
associated zero-curvature equations. At last in the fourth section we consider integrable hierar- 
chies associated with loop algebras in different gradings and different i2-operators. 

2 Commutative subalgebras and classical R operator 

Let g be a Lie algebra (finite or infinite-dimensional) with a Lie bracket [ , ]. Let -R : g ^ g 
be some linear operator on g. The operator R is called the classical i?-operator if it satisfies 
modified Yang-Baxter equation: 



We will use, in addition to the operator R, the following operators: R± = i? it 1. As it is 
known [21 [12] the maps R± define Lie subalgebras qr^ C g: Qr_^ = Imi?-|-. It is easy to see from 
their definition gjj^ + qr = g, but, in general, this sum is not a direct sum of vector spaces. 

Let g* be the dual space to g and ( , ) : g* x g ^ C be a pairing between g* and g. 
Let {Xi} be a basis in the Lie algebra g, {X*} be a basis in the dual space g*: (X*, Xi) = 6ij, 
L = LiX^ € g* be the generic element of g*, Lj be the coordinate functions on g*. 




RmX),Y] + [X,R{Y)]) - [R{X),R{Y)] = [X,Y] 



vx,y eg. 



(2) 



Let us consider the standard Lie-Poisson bracket between Fi,F2 G C°°(g*): 



{Fi(L),F2(L)} = (L,[VFi,VF2]), 



where 




is a so-called algebra- valued gradient of Ff^. The summations hereafter are implied over all basic 
elements of g. 
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Let R* be the operator dual to R, acting in the space q*: {R*{Y),X) = {Y,R{X)), G 5*, 
X G g. Let I^{q*) be the ring of invariants of the coadjoint representation of q. 
We will consider the functions (L) on g* defined by the following formulas: 

I^^{L) ^ U{R* ±1){L)) ^ UR^L)). 

The following theorem holds true |13j : 

Theorem 1. Functions {/^"""(L)} and {ll^'{L)}, where Ik{L),Ii{L) G generate an 

Abelian subalgebra in C°°{q*) with respect to the standard Lie-Poisson brackets { , } on q* : 

(i) {I^^L),lf^L)} = 0, (n) {/f-(L),/f-(L)} = 0, (m) {/f (L), if " (L)} = 0. 

Remark 1. Note that the commutative subalgebras constructed in this theorem differ from 
the commutative subalgebras constructed using standard i?-matrix scheme [12]. Indeed, our 
theorem states commutativity of functions {/^"''(L)} and {/^ ~ (L)} with respect to the initial 
Lie-Poisson bracket { , } on g. The standard i?-matrix scheme states commutativity of the 
functions {/^(L)} with respect to the so-called i2-bracket { , where: 

{Fi{L),F2{L)}r = (L,[i?(VFi),VF2] + [VFi,/2(VF2)]). 

Theorem [T] provides us a large Abelian subalgebra in the space {C°°{q*), { , }). We will con- 
sider the following two examples of the i?-operators and the corresponding Abelian subalgebras. 

Example 1. Let us consider the case of the Lie algebras g with the so-called "Kostant-Adler- 
Symes" (KAS) decomposition into a direct sum of the two vector subspaces g±: 

5 = 0+ + 0-, 

where subspaces g± are closed Lie subalgebras. Let P± be the projection operators onto the 
subalgebras q± respectively. Then it is known [12] that in this case it is possible to define the 
so-called Kostant-Adler-Symes iZ-matrix: 

i? = P+ -P_. 

It is easy to see that R+ = 1 + R = 2P+, ii_ = R — 1 = — 2P„, are proportional to the projection 
operators onto the subalgebras g±. It follows that = q± and gR_^ <^Qr_ = 0. 
The Poisson commuting functions (L) acquire the following simple form: 

I^^{L) = I^{L) = h{L±), where L± = P^L, 
i.e. I^{L) are restrictions of the coadjoint invariants onto the dual spaces g^. 
Example 2. Let us consider the case of Lie algebras g with the "triangular" decomposition: 

5 = 0+ + 00 + 0-, 

where the sum is a direct sum of vector spaces, g± and go are closed subalgebras, and g± are 
go-modules. As it is known [4] (see also [13] for the detailed proof), if Rq is a solution of the 
modified Yang-Baxter equation dJj) on go then 

R = P+ + Ro-P^ (3) 

is a solution of the modified Yang-Baxter ^ equation on g if P± are the projection operators 
onto the subalgebras g± . 
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In the case when Rq = HcIq (which are obviously the solutions of the equation ^ on go) we 
obtain that /^-matrix ^ passes to the standard Kostant-Adler-Symes i?-matrix. Nevertheless 
in the considered "triangular cases" there are other possibilities. For example, if a Lie subalgebra 
00 is Abelian then ([3]) is a solution of ([2]) for any tensor Rq on qq. 

In the case of the i?-matrix ([3]) we have: R+ = 2P+ + (Pq + ^o), R~ = -(2P_ + (Pq - Rq)), 
9R± = 0± + I™ {Ro)±, where {Ro)± = {Ro ± Pq) are the i?±-operators on go and Qn+ H qr_ = 
Im (i?o)+ n Im (i?o)-- The Poisson-commutative functions (L) acquire the following form: 



I^(L)./f-(L) = 4(Li + ii±^(Lo)) 



where L± = P^L, Lq = PqL. We will use such functions constructing different generalizations 
of the Thirring model. 

The Theorem [T] gives us a set of mutually commutative functions on g* with respect to the 
brackets { , } that can be used as an algebra of the integrals of some Hamiltonian system 
on g*. For the Hamiltonian function one may chose one of the functions I^^{L) or their linear 
combination. Let us consider the corresponding Hamiltonian equation: 

^ = {L.,lf(L)}. (4) 
The following proposition is true |13j : 

Proposition 1. The Hamiltonian equations of motion (jj]) can he written in the Euler- Arnold 
(generalized Lax) form: 

^ = ad* ±L, (5) 
where V^ = VI^^iL). 

Remark 2. In the case when g* can be identified with g and a coadjoint representation can 
be identified with an adjoint one, equation ([5]) may be written in the usual Lax (commutator) 
form. 

In the present paper we will not consider finite-dimensional Hamiltonian systems that could 
be obtained in the framework of our construction but will use equation ([5]) in order to generate 
hierarchies of soliton equations in 1 + 1 dimensions. 

3 Integrable hierarchies and classical i?-operators 

Let us remind one of the main Lie algebraic approaches to the theory of soliton equations [10] . 
It is based on the zero-curvature condition and its interpretation as a consistency condition of 
two commuting Lax flows. The following Proposition is true: 

Proposition 2. Let Hi,H2 G C°°(g*) be two Poisson commuting functions: {Hi,H2} = 0, 
where { , } is a standard Lie-Poisson brackets on g*. Then their q- valued gradients satisfy the 
"modified" zero- curvature equation: 

where I is a Casimir function and ti are parameters along the trajectories of Hamiltonian equa- 
tions corresponding to the Hamiltonians Hi and k is an arbitrary constant. 
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For the theory of sohton equations one needs the usual zero-curvature condition (case k = 
in the above proposition) and the infinite set of the commuting Hamiltonians generating the cor- 
responding U-V pairs. This can be achieved requiring that g is infinite-dimensional of a special 
type. In more detail, the following theorem is true: 

Theorem 2. Let q be an infinite- dimensional Lie algebra of Q-valued function of the one complex 
variable u and a Lie algebra g be semisimple. Let R be a classical R operator on g, L{u) be the 
generic element of the dual space q* , Lk{L{u)) be Casimir functions on g* such that functions 

D , _____ D , 

{L{u)) are finite polynomials on g*. Then the Q-valued functions V/^ {L{u)) satisfy the 
zero- curvature equations: 

Proof. As it follows from Proposition [2] and Theorem [1] the algebra-valued gradients 
V/f±(L(u)), V/f±(L(n)) and V/f±(L('u)), V/f^(L(n)) satisfy the "modified" zero-curvature 
equations On the other hand, as it is not difficult to show, for the case of the Lie algebras g 
described in the theorem the algebra- valued gradients of the Casimir functions are proportional 
to powers of the generic element of the dual space L{u), i.e. to the formal power series. On the 
other hand, due to the condition that all {L{u)) are finite polynomials, their algebra- valued 
gradients are finite linear combinations of the basic elements of the Lie algebra g. Hence the 
corresponding modified zero-curvature equations are satisfied if and only if the corresponding 
coefficients k in these equations are equal to zero, i.e. when they are reduced to the usual 
zero-curvature conditions. This proves the theorem. ■ 

Remark 3. Note that equations ([7]), ^ define three types of integrable hierarchies: two "small" 
hierarchies associated with Lie subalgebras gR_,_ defined by equations d?]) and one "large" hierar- 
chy associated with the whole Lie algebra g that include both types of equations ([7]) and ([8]). 
Equations ([8|) have an interpretation of the "negative flows" of the integrable hierarchy asso- 
ciated with Qrj^. In the case g/j^ ~ g/j_ the corresponding "small" hierarchies are equivalent 
and the large hierarchy associated with g may be called a "double" of the hierarchy associated 
with qr^ . 

In the next subsection we will consider in detail a simple example of the above theorem when 
g is a loop algebra and i?-operator is not of Kostant-Adler-Symes type but a triangular one. We 
will be interested in the "large" hierarchy associated with g and, in more detail, in the simplest 
equation of this hierarchy which will coincide with the generalization of the Thirring equation. 



4 Loop algebras and generalized Thirring equations 

4.1 Loop algebras and classical i?-operators 

In this subsection we remind several important facts from the theory of loop algebras |8]. 

Let g be semisimple (reductive) Lie algebra. Let g = ^ gj be Zp = Z/pZ grading of g, i.e.: 

3=0 

[BIjSj] C g^qTj where j denotes the class of equivalence of the elements j ^ 1> modp Z. It is 
known that the Zp-grading of g may be defined with the help of some automorphism a of the 
order p, such that o"(gj) = e^'^^^/^Qj and gjy is the algebra of a-invariants: 0"(gQ) = gg. 
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Let = 0® Pol(ti, u ^) be a loop algebra. Let us consider the following subspace in 0: 
0" = 00j0n^'. 

It is known [5] that this subspace is a closed Lie subalgebra and if we extend the automorphism a 
to the map a of the whole algebra 0, defining its action on the space (8) Pol(u, in the 
standard way [^r a{X (g) u^) = (t{X) (g) g-^Tifc/pyA;^ i^a.eTL the subalgebra J'^ can be defined as the 
subalgebra of a-invariants in 0: 

= {X® p{u) G 0| a{X (g) p{u)) =X® p{u)]. 

We will call the algebra 0'^ the loop subalgebra "twisted" with the help of a. The basis in 0'^ 
consists of algebra- valued functions {X^ = XaU^}, where Xa € 0j. Let us define the pairing 
between 0^^ and (0°^)* in the standard way: 

{X,Y) = Tesu"\X{u),Y{u)), 

u=0 

where X £ 0'^, Y € (0'^)* and ( , ) is a bilinear, invariant, nondegenerate form on 0. It is easy 
to see that with respect to such a pairing the dual space (0^^)* may be identified with the Lie 
algebra 0*^ itself. The dual basis in (0^^)* has the form: {Yi = X~^'°'u~^}, where X~^'°' is a dual 
basis in the space 03j. 

The Lie algebra 0*^ possesses KAS decomposition 0°" = 0"^^^ + 0*^" flT], where 

= 0J ® ' 5""" = 9j ® 

It defines in a natural way the Kostant-Adler-Symes ii-operator: 

R = P+ -P-, 

where are projection operators onto Lie algebra 0^^^. 

The twisted loop algebra 0^^ also possesses "triangular" decomposition: 

r = 0^+00+0^, 

where g"- = q'L, 0'"+ =0^ + 00. 

It defines in a natural way the triangular i?-operator 

R = P+ + Rq — P_ , 

where P± are the projection operators onto Lie algebra 0o-± and Rq is an /^-operator on 0o. 

The Lie algebra 0o in this decomposition is a reductive Lie subalgebra of the Lie algebra 0. 
Due to the fact that a lot of solutions of the modified classical Yang-Baxter equations on 
reductive Lie algebras are known, one can construct explicitly /^-operators Rq on 0o, the Lie 
subalgebras qr^ = 0± + Iin(i2o)± and the Poisson-commutative functions constructed in the 
Theorem [TJ 

4.2 Generalized Thirring models 

As we have seen in the previous subsection, each of the gradings of the loop algebras, corre- 
sponding to the different automorphisms a yields its own triangular decomposition. Let us 
consider the simplest equations of integrable hierarchies, corresponding to different triangular 
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decompositions in more detail. The generic elements of the dual space to the Lie algebras qfi.± 
are written as follows: 

/ dimgg \ ±00 dimgj 



\ a=l j j=itl 0=1 i=±l 

where L^~^^ G and in order to simplify the notations we put r± = -^i^^. 

Let ( , ) be an invariant nondegenerated form on the underlying semisimple (reductive) Lie 
algebra q. Then hiL) = \{L,L) is a second order Casimir function on q* . The corresponding 
generating function of the second order integrals has the form: 



lt{L{u))=h 



/ ±CXD ±00 \ 

(L±(n)) = - + ^ L(^■)n-^ri(LW) + ^ 

\ J=±i i=±i / 



±00 



where p is an order of a. The simplest of these integrals are: 



The algebra-valued gradients of functions if read as follows: 

p-i 

Vlf{u) = n±P4(LW) +Y^u^(p-^)l^±^) + r±(L(±P)). (9) 

i=i 

They coincide with the U-V pair of the generalized Thirring model. The corresponding zero- 
curvature condition: 

avi^ _ awgLM) , [v,j.(„„),. w-,L,„,), . (10) 

yields the following system of differential equations in partial derivatives: 

-5,„r;(L(0)) = [r;(L(0)),r„(L("^'))], (11a) 
a,^r_(L(°)) = [r+(L(P)),rl(LW)], (lib) 

= [L(?'-^'),r_(L(-P))] + E + [rUL^^^), L^"% (11c) 

p-i 

5^^L(^'-p) = [r+(L(p)),L(^-P))]+ ^ 

+ [L(^-),rl(LW)],jGl,p-l, (lid) 

p~i 

9,^r„(L(-P)) -a,_r+(L(P)) = [r;(L(o)), rl(LW)] + E[^^''"'^ i^'"''^] 

1=1 

+ [r+(L(P)),r_(L(-P))]. (lie) 

These equations are integrable generalization of Thirring equations corresponding to Lie algeb- 
ra 3, its Zp-grading defined with the help of the automorphism a of the order p and the classical 
i?-operator Rq on Qq. We will call this system of equations the non-Abelian generalized Thirring 
system. In order to recognize in this complicated system of hyperbolic equations generalization 
of Thirring equations we will consider several examples. 
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4.2.1 Case of the second-order automorphism 

Let us consider the case when automorphism a is involutive, i.e. p = 2. Then g- = g— and Lax 
matrices L±(u) are the following: 

L±(n) = ri(LW) + u^^L^^^^ + u^^L^^^^ + • • • , 

and the U-V pair ([9]) acquire more simple form: 

Vlt\u) = n±24(LW) + n±^L(±i) + r±(L(±2))_ (12) 

The corresponding zero-curvature condition yields the following system of differential equations 
in partial derivatives: 

a._r;(LW) = [r_(L(-2)),r;(L(0))], (13a) 

a,^r:l(L(0)) = [r+(L(2)),r:l(LW)], (13b) 

a,_L« = [r_(L(-2)),L(i)] + [L(-i),r;(LW)], (13c) 

a,^L(-i) = [r+(L(2)),L(-i)] + [L(i),r:l(L(0))], (13d) 
5.,r„(L(-2)) _ 5^_,^(^(2)) ^ [,;(l(0)),,1(lW)] + 

+ [r+(L(2)),r_(L(-2))]. (13e) 



The system of equations (I13p is still sufficiently complicated. In order to recognize in this 
system the usual Thirring equations we have to consider the case q = sl{2). 

Example 3. Let q = sl{2) and a be the Cartan involution, i.e. sI{2)q = f) = diag(a,— a) and 
sl{2)i consists of the matrices with zeros on the diagonal. The Lax matrices L±{u) have the 
following form: 



-a(0) J^"" [ 7(±i) j V -a(±2) 

Due to the fact that sI{2)q is Abelian the one-dimensional linear maps r± are written as follows: 
r'^{L^^^) = k±L^^^ where k± are some constants, such that k+ + k- = 1. 
The simplest Hamiltonians obtained in the framework of our scheme are: 

If = A:±(a(0))2, If = 2A:±a(0)a(±2) + 

As it follows from the Theorem [1] and, as it is also easy to verify by the direct calculations, these 
Hamiltonians commute with respect to the standard Lie-Poisson brackets on sl{2) : 

|Q,{2fc) ^ p{2l+l) I ^ ^{2{fc+0+l) ^ |Q,{2fc) ^ ^{2/+l) I ^ _^(2{fc+0+l) ^ 

{^(2fc+l)^^{2/+l)| ^ 2a(2(fc+0+2)^ 

{a(2fc),a{20} = |^(2fc+l)^^(2m)| ^ |^(2A:+1)^^(2«+1)| ^ q_ 

As a consequence, they produce the correct U-V pair for zero-curvature equations: 
±2 _ , a^^'^ \ ±1 / A , ±2 



where C/ = VI f, V = VI f. 

Due to the fact that if are constants of motion we have that a^^^ is also a constant of 
motion and equations (|13ap . ()13bp are satisfied automatically. The integrals if are constants 
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of motion too, that is why a^^^^ are expressed via 7^^^^ Hence, equation (jl3ep is not 

independent and follows from equations (|13c|) . (|13dp . These equations are independent. The last 
summand in the equations (jl3cp . ()13dp gives the linear "massive term" in the Thirring equation. 
The first summand gives the cubic non-linearity. 

Let us write the equations (jl3cp . ()13dp in the case at hand in more detail taking into account 
the explicit form of the matrices L^'^^ and linear operators r±: 

a,._7^'^ = 2fc+a(°)7^-'^ - 2fe_a(-2)7(i), 
5,^7(-i) = 2fc_aW7(^) - 2fc+a(2)7{-i), 

where we have put I^^ = and, hence, a^^^^ = ~ 2fc^a(o) (Z^^^^-^T^^^^)- 

The system of equations (|14p admits a reduction a^^^ = ic, 7^='=^) = —f3^^^^ = — V'ii, which 
corresponds to the restriction onto the real subalgebra su{2) of the complex Lie algebra sl{2). 

After such a reduction the system of equations (jl4p is simplified and acquires the form: 

In the case = k- = ^ these equations are usual Thirring equations with a mass m = c^. 
4.2.2 Matrix generalization of Thirring equation 

Let us return to the generalized Thirring model in the case of the higher rank Lie algebra g, 
its automorphism of the second order and corresponding Z2-grading of g: = 0o + fli- We are 
interested in the cases that will be maximally close to the case of the ordinary Thirring equation. 
In particular, we wish to have r^{L^^^) that enter into our U-V pair (jl2p near the second order 
of spectral parameter to be constant along all time flows. As it follows from the equations (jl3ap . 
()13bp , this is not true for the general r-matrix Rq on gg . Fortunately there are very special cases 
when it is indeed so. The following proposition holds true: 

Proposition 3. Let Qq admit the decomposition into direct sum of two reductive subalgebras: 
00 = Bq' © 0g ■ = + ^-^ be the corresponding decomposition of the element of the 

dual space. Let C(0g^) a center of the subalgebra g^, K{L^^) be a part of L^^ dual to the 
center of the subalgebra 0^. Then Rq = Pq — is the R operator on Qq, r^{L^^^) = L^_^\ 

r^{L^^^) = L?"* are constant along all time flows and the reduction L^^^ = Kj^ + where 
K± is a constant element of C(0g^)*; ^-^ consistent with all equations of the hierarchy d?]), ^ 
corresponding to the loop algebra and triangular R-operator on 0'^ with the described above 
KAS R-operator Rq on Qq. 

Proof. Let us at first note that in the case under consideration r^ = r±. Let us show that 
r^(L(^)) = r±{L^^^) are constant along all time flows generated by the second order Hamilto- 
nians L^"^^. It is easy to show that the corresponding algebra- valued gradients have in this case 
the following form: 

where P± are the projection operators onto the subalgebras 0^ in the triangular decomposition 
of the Lie algebra 0°^ and we took into account that r+-|-r_ = 1. Let us substitute this expression 
into the Lax equation: 



dt^ 



- = [WL}'\L{u)], (15) 

2k 
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and take into account that the case under consideration corresponds to the Kostant-Adler-Symes 
decomposition and, hence, we have correctly defined the decomposition L{u) = L^{u) +L„(n), 
and each of two infinite-component Lax equation (jlSp is correctly restricted to each of the 
subspaces L±{u): 



Making projection onto the subalgebra Qq in these equations we obtain the equations: 
d^:^ = -[r^(L(±2'=)),r±(L(0))], = [r^(L(^2'=)), r^lL^)]. 

"-hk "''2/t 

Due to the fact that our r-matrix is of Kostant-Adler-Symes type, we have [r^{X), r±{Y)] = 0, 
y X,Y G gg and, hence, r+(L^'^^) = r_(L^*'^) = are constant along all time flows 

generated by In analogous way it is shown that L^^^ = + are constant with 

respect to the time flows generated by the higher order integrals I^^^- Hence, in this case 
components of L^^^ belong to the algebra of the integrals of motion of our infinite-component 
Hamiltonian or Lax system. This algebra of integrals is, generally speaking, non-commutative 
i.e. {la\l^^^} = where c^^ are structure constants of the Lie algebra gg- That is why 

in order to have a correct and consistent reduction with respect to these integrals we have to 
restrict the dynamics to the surface of zero level of the integrals belonging to ([flOiflo])*- Other 
part of (so)*, namely the one belonging to = (C(0o))* = (C(0o ))* + (C(0o ))*, may 

be put to be equal to a constant, i.e. correct reduction is: L^^'^ = + J^_, where K± is 
a constant element of C(0g^)*- ^ 

Hence, in this case we have the following form of simplest U-V pair ()12p for the zero-curvature 
condition: 

Vlt\u) = u^^Ki + + (16) 

where L^^-* = P^{L^^'^'^) e g^. The corresponding equations (I13ap . (|13b|) are satisfied auto- 
matically and the rest of equations of the system (|13p are: 

= + [l(-'\k+], (17a) 

a,.^L(-^) = [Lf ,L(-i)] + (17b) 

= [lW,L(-^)], (17c) 

where we have again used that in our case [r+(X), r_ (K)] = 0, yX,Y € %. 

The equations (jl7ap . (|17bp are the simplest possible generalizations of the Thirring equations. 
The first term in the right-hand-side of this equation is an analog of the cubic non-linearity of 
Thirring equation. The second term is an analog of linear "massive" term of the Thirring 
equations. In all the cases one can express L^'^^ via L^^^\ and equation (jl7cp will follow from 
equations (I17ap . (jl7bp . In order to show this we will consider the following example. 

Example 4. Let us consider the case g = gl{n), with the following Z2-grading: g = gl{n)Q + 
gl{n)i, where gl{n)Q = gl{p) + gl{q) and gl{n)i = C^^'^, i.e. 

5/(n)o=(o ]). 5Kn)i=(° I 
where a € gl{p), 6 G gl{q), (3 € Mat(p, g), 7 € Mat(g,p). 
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In this case Qq = gl{p), flg = gl{q), C(flg ) = k+lp, ((flg ) = ^-^Q- '^^^ corresponding U-V 
pair (|16|) has the form: 

The corresponding zero-curvature condition yields in this case the following equations: 

d^J+ = -0+6. + kj.), a,_7+ = (L7+ + A:+7_), (18a) 
d,^i3- = {a+i3- + /._/?+), 9.+7_ = -(7-a+ + A:-7+), (18b) 
= (7+/3- - 7-/3+), 9,_a+ = (/3_7+ - 7-/3+)- (18c) 

By direct verification it is easy to show that the substitution of variables: 

1 . 1 . 

(5_ = -— 7_/3_, a+ = -— /3+7+ 

solves equation (jl8cp and yields the following non-linear differential equations: 

d.J+ = ■^/5+(7-/3-) - k+p., a,_7+ = -■^(7-/3-)7+ + Kl-, 
dx+P- = -■^(/3+7+)/3- + k-(3+, dx^^- = ■^7_(/3+7+) - A;_7+. 

These equations are matrix generalization of the complex Thirring system |15j . They admit 
several reductions (restrictions to the different real form of the algebra gl(n)). For the case of 
an u{n) reduction we have 

7± = -PI. = ^± = ^^±, G M 

and we obtain the following non-linear matrix equations in partial derivatives: 

1 1 

= — ^'+(^'L^'_) + «:+^'_, ^_ = (^'+1'\.)^'_ + 

«;__ K+ 

These equations are a matrix generalization of the massive Thirring equations ( ^± E Mat(p, g)). 
In the case ofp = n — 1, q = 1 or p = 1, q = n — 1 we obtain a vector generalization of Thirring 
equations. In the special case n = 2,p = q = l and k_ = k+ we recover the usual scalar massive 
Thirring equations with mass m = (k+)^. 

4.2.3 Case of Coxeter automorphism 

Let fj be a Coxeter automorphism. In this case p = h, where /i is a Coxeter number of g 
and algebra is Abelian. That is why all tensors Rq are solutions of the mYBE on and 
maps r± are arbitrary (modulo the constraint r+ -|-r_ = 1). Moreover, in this case the following 
proposition holds: 

Proposition 4. Let a be a Coxeter automorphism of g and maps r± on g be nondegenerated. 
In this case L^^^ is constant along all time flows and components of L^^^^ are expressed as 
polynomials of components of L^^^\ . . . , L^^^'^'^^^ 
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Proof. In the considered case we have that dirngg = rankg. On the other hand there exist 
r = rankg independent Casimir functions IkiL), k € l,/i, on q*. Hence there exist 2rankg 
integrals of the following form: 

=4(ri(LW)). 

They are constant along all time flows generated by all other integrals I^K That is why we 
may put them to be equal to constants. On the other hand, it is not difficult to see that 
all r independent components of L^^'^ can be functionally expressed via I^^. Hence, they are 
constants too. At last, we have 2 rankg integrals I^^, which are constant along all time flows 
and we may put them to be equal to constants I^^ = const^'^. It easy to see that I^'^ are linear 
in L^^'^) and, on the surface of level of I^^, all components of L^^'^) are expressed polynomially 
via components of L^^^\ . . . , L^^^^~^^^ if the maps r± are non-degenerate. ■ 

This proposition has the following important corollary: 

Corollary 1. The number of independent fields in the generalized Thirring equation (llOp . cor- 

h-l 

responding to Coxeter automorphism, is equal to 2 ^ dimg^- = 2(dimg — rankg). 

i=i 

Let us explicitly consider Thirring-type equations (jlip in the case of the Coxeter automor- 
phisms. In this case gjj ~ g^^ is Abelian and equations (jllap . (jllbp become trivial. Moreover, due 
to the fact that L^^'^-* is expressed via L^^^^ where j < h, equation (jlldp becomes a consequence 
of equations (jllcp , (jlldp . In the resulting system of equations ([TT]) is simplified to the following 
system: 

i=j+i 

d^^L^^-^^ = [r+(L('^)),L(j'-'^))] + + [L(^'),r* (L(°))], (19) 

i=j+i 

where j G {l,/i - 1}, is a constant matrix and L^^^^ is polynomial in L^^^\ We call 
this system of equations the Abelian generalized Thirring equations. The last summand in 
equations (jl9p is an analog of the "massive term" in the Thirring equation. The first summand 
is an analog of the cubic non-linearity in the Thirring equation, but with the growth of the 
rank of g the degree of this term is also growing. The other terms are of the second order in 
dynamical variables. They are absent in the ordinary Thirring system corresponding to the case 
g = sl{2). Let us explicitly consider the simplest example, which already possesses all features 
of the generalized Thirring equation: 

Example 5. Let g = gl{3). Its Za-grading corresponding to the Coxeter automorphism has the 
following form: 

/qiOO\ /0/3iO\ /00 73 \ 

5^(3)5 = a2 , gmi = /?2 , g^h = 7i , 
\ as J \ (33 J VO72O/ 

5^(3)3 = 5^(3)0, 5^(3)31 = 9^(3)2, 5^(3)32 = am-i- 

The Lax operators belonging to the dual spaces to gl{3)j^^ are: 

L±(n) = ri(LW) + L^^'^u^' + L^^^^u^^ + L(±3)^t3 + . . . ^ 
where L(o),L(±3) g g/(3)g, g gl{3)-2, L^^\L^~^^ G gl{3)i. 
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In order to simplify the form of the resulting soliton equations we will use the following 
notations: 
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The generating functions of the Poisson-commuting integrals of the corresponding integrable 
system are: 



/fc(L±(u)) = itr(L±(M))^ = Ik'^u^""- 



k 

m=0 

In particular, we have the following integrals: 
If = ^tr(4(LW))^ fceM. 

From the fact that they are fixed along all the time flows we obtain that all ctj are constants of 
motion. We also have the following integrals: 

If = trL^^'\ If = tr (4(lW)L(±3)) + tr 

If = - tr (L(±i))3 + tr (4(lW)(L(±i)l(±2) + L^±^)l^±^))) + tr {rUL^^^ f L^^'^) . 
3 

These integrals permit us to express components of L^^^^ via components of L^^^^ and L^^^^ in 
the case of the non-degenerated maps r±. Let us do this explicitly. For the sake of simplicity 
we will put that r'j_{L^'^^) = k±L^^\ where k± are some constants and /c+ + /c_ = 1. In this case 
we will have the following explicit form of the above integrals: 



i=l i=l i=l 

3 3 



1=1 1=1 

where we have implied in the last summation that (3^ = (3^ , 7^ = 7^ . 
The U-V pair corresponding to the Hamiltonians if have the form: 

Vlf{u) = u±34(L(0)) + n±2L(±i) + ^,±1l(±2) + r±(L(±3)). 
The corresponding zero-curvature equation reads as follows: 

5.,L(-^) = [r+(L(=^)),L(-i)] + [L(2),rl(LW)], 

d,^L^~^^ = [r+(L(3)),L(-2)] + [L(i),rl(L(0))] + [L^^),L^-% 
In the component form we have the following equations: 

dx.Pt = k-i^i - ^i+i)(^t + k+{ci+i - Ci)7j", 
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d^^lt = KiSr^, - + k-ic^ - Q+i)/?r + _ j+p-), (20) 

where q = Ofj, indices i, j, k constitute a cyclic permutation of the indices 1, 2, 3, and it is 
implied that 6^^^ = 6^, C3+1 = ci. 

Taking into account that all the integrals are constants of motion and putting their values 
to be equal to zero we can explicitly express variables 5f via Oj = q, P^, 

1 



k'j^ici Cj){ci Cfc) 

X (k±{c,+Ck)Y,P^^t-k±Y.ci{P^l^+Pti-fti)-PtPtPt] , (21) 

V 1=1 1=1 / 

where indices j and k are complementary to the index i in the set {1,2,3}. Due to the fact 
that Ci are constants of motion this expression is polynomial of the third order in dynamical 
variables. 

At last, substituting (j2ip into the equations (I20p we obtain a system of nonlinear differential 
equations for the dynamical variables and j^. These equations are a 5/(3) generalization 
of "complex" Thirring equations. Unfortunately, neither these equations nor more compli- 
cated gl{n) complex Thirring equations do not admit reductions to the real forms u{3) or u{n) 
respectively. 



4.3 The case of Kostant— Adler-Symes i2-operators 

Let us now consider integrable hierarchies corresponding to the case r_|_ = 1, r_ = (or r„ = 0, 
r_ = 1), i.e. corresponding to the Kostant -Adler-Symes decomposition of t loop algebras. 

As it was reminded in Subsection l4.lt each of the gradings of the loop algebras corresponding 
to the different automorphisms a of order p provides its own Kostant-Adler-Symes decompo- 
sition: q'^ = + S*^", and, hence, provides commuting Hamiltonian flows and hierarchies of 
integrable equations. 

Let us consider the simplest equations of integrable hierarchies corresponding to them. We 
have the following generic elements of the dual spaces (0°^^)*: 



i=i i=i 



where ^ 0fc) ^-rid the following generating functions of the commutative integrals: 

00 00 

4(L+(n)) = iru-^, ii{L-{u)) = Y.iru\ 

m=0 n=l 

where Ik{L) are Casimir functions of g. 

We will be interested in the following Hamiltonians and their Hamiltonian flows: 

If = (lW,lW) + i'^(L(^),L(P-^)), I^P = /,(L("i)). 
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The algebra- valued gradients of functions , Ip^ read as follows: 

p-i 

i=i 

where L«= E ffeX"! G g- 
i=i 

The corresponding zero-curvature condition yields the following system of equations: 

5,._L(°)=0, (22) 
a,_L('=) = [Z«,L(^-i)], fcG{l,p}, (23) 
a,^LW = [L(^'),LW]. (24) 

Remark 4. The gradient V/|(n) is an analog of the U operator of the generalized Thirring 
hierarchy. The other gradient VIp^{u) is an analog of the V operator of the non-Abelian Toda 
equation. Hence the corresponding integrable equations may be considered as an intermediate 
case between generalized Thirring and (non-Abelian) Toda equation. 

4.3.1 Case of second order automorphism 

Let us consider the case of the automorphism of the second order o"^ = 1 (j) = 2), the correspon- 
ding Hamiltonians, their matrix gradients and zero-curvature conditions. We will use commuting 
second order Hamiltonians of the following form: 

I| = (lW,L(2)) + i(LW,LW), = i(L(-i),L(-i)). 

The algebra- valued gradients of functions /I) read as follows: 

Vl|(n) = 7.2^(0) +nL(i) + L(2), V/-P(n) = u-^L^-^\ 

The corresponding zero-curvature condition yields the following simple system of differential 
equations of hyperbolic type: 

a.,L("^) = [L(2),L("^)], (25) 

where g go, ^(±1) G gj. 

Let us consider the following example of these equations: 

Example 6. Let g = sl{2) and cr be a Cartan involution, i.e. s/(2)q = diag(Q,— a) and sZ(2)j 
consists of the matrices with zeros on the diagonal. The simplest Hamiltonians obtained in the 
framework of Kostant-Adler-Symes scheme are: 

/o = (a{0))2, j2 ^ 2aWa(2) +^(1)^(1)^ j2 ^ ^(-i)^(-i). 
They produce the following U-V pair for zero-curvature equations: 



V/l = ( „ \+u[ '\ \+u 



a(2) 




















\ 2 ( 



2 - I -a(2) ; ^ " V 7^^) ; ^ " I -aW 
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Due to the fact that I2 are constants of motion we have that a^'^^ is also a constant of motion. 
Moreover, using the fact that /| is a constant of motion one can express a^'^^ via 7(1). 
Hence we obtain from the equations (j25p the following independent equations for the variab- 
les /?(±i),7(±i): 

5^.7(1) = 2a(°)7(-i), (26a) 
= -2a(°) (26b) 
5^+7(-i) = -2a(2)^(-i)^ (26c) 
5,^/3(-i) = 2a(2)/3(-i). (26d) 

From the equation ()26bp it follows that Z^^"-*^) = — '^^^(o) • Substituting this into equation ()26dj) 

and taking into account that a^*^^ is a constant along all time flows we obtain the following 
equation: 

'9^+._/3^'^ = 2a(2)a,_/?«. (27) 

Taking into account that a^^^ = — ^^io) where we have put that /| = 0, and making 

the reduction to the Lie algebra su{2): a^^^ = ic, 7^"^^ = —(3^^^ = — we obtain the following 
integrable equation in partial derivatives: 



5'+x_^ + -IV'l'5,_V' = o. 



This equation is (in a some sense) intermediate between Thirring and sine-Gordon equations. 



4.3.2 Case of Coxeter automorphism (principal gradation) 

Let us consider again the case of the principal gradation. In this case p = h and a subalgebra gg 
is Abelian. In the same way as it was done in the case of the "generalized Abelian Thirring 
models" it is possible to show that L^^^ is constant along all time flows and components of L^^^ 
are expressed polynomially via the components of L^*^) , k < h. Let us assume that constants of 
motion L^''^ are such that the operator adj;^{o) is nondegenerate. In such a case we may solve the 
first of the equations (j23p in the following way: 

L« = -ad-(^o)(5._L(')). 

Substituting this expression into equation (j24|) and taking into account commutativity of gg and, 
hence, operators ad^^gj and adj;^{h) , we finally obtain the following matrix differential equation 
in partial derivatives: 

5.,._L« = [lW(L«,...,l('^-i)),5._L«], 
where L^''\ A; G {2, /i — 1} satisfy the following set of ordinary differential equations: 

a,._LW = [L('^-i),ad-(^,(a,_L«)]. 

Let us consider the following example. 

Example 7. Let g = gl{3), h = 3. In this case we have the following differential equations: 

5.,._L« = [L(3)(L«,L(2)),9._lW], (28) 
5,_L(2) = [L«,ad7(^„)(a,_L«)]. 
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Let lW, lW, L(2), are parametrized as in the Example O i.e. 



ai 
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1 , L(^^ = 
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In such coordinates equation (I28p acquires the fohowing form: 

Cx+A = ('5m(A7) - <5^(/3,7))5x_ft, i e {1,3}. (29) 
where, hke in the Example O 7) are expressed via (3j, 7^ and constants ai = cf. 

Si{f3, 7) = 7 7 I (9 + cfc) > J A7« - > : ciiPni + A-17Z-1) - /3i/?2/33 1 • (30) 



{a - cj){ci - cfc) 



'33 \ 

{cj + Cfc) A7« - Yl ^^^^i^i + A-i7z-i) - 



This equation is an exact analog of the equation (|27p . Unfortunately in this case there is no 
su{3) reduction and variables ji are not conjugated to Pi but satisfy the following differential 
equations: 

dx-li = eijkPj{ck - Ck+i)~^dx_l3k, i e {1,3}. (31) 

Equations (|29|) ~ (j31|) are intermediate between the Abelian 5/(3)-Toda equations and the genera- 
lized Abelian (7/(3)-Thirring equations. 
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